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ABSTRACT 

NGC 4194 is a post-merger starburst known as The Medusa for its striking tidal features. We present here a detailed 
study of the structure and kinematics of ionized gas in the central 0.65 kpc of the Medusa. The data include 
radio continuum maps with resolution up to O'/ 18 (35 pc) and a 12.8 /zm [Nen] data cube with spectral resolution 
~4 km s -1 : the first high-resolution , extinction-free observations of this remarkable object. The ionized gas has 
the kinematic signature of a core in solid-body rotation. The starburst has formed a complex of bright compact 
H ii regions, probably excited by deeply embedded super star clusters, but none of these sources is a convincing 
candidate for a Galactic nucleus. The nuclei of the merger partners that created the Medusa have not yet been 
identified. 

Key words: galaxies: individual (NGC 4194) - galaxies: kinematics and dynamics - galaxies: starburst - galaxies: 
star formation 
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1. INTRODUCTION 

The most luminous starbursts appear to be caused by major 
mergers of galaxies (e.g., Sanders et al. 1988) and are rare 
because major mergers are infrequent. Minor mergers create 
LIRGS; slightly less luminous but much more common, and 
important to galaxy evolution due to their sheer numbers 
(Kartaltepe et al. 2010). Minor mergers, like major mergers, 
appear to instigate intense star formation in compact regions 
and create concentrated sources. But minor mergers evolve 
differently from major mergers, so the characteristics of star 
formation are likely to be different. How does star formation 
develop in a minor merger? How will these starbursts evolve 
and how do they affect their surroundings as they contribute to 
the newly formed galaxy? To answer those questions we need 
to understand the structure and kinematics of the gas ionized 
by the embedded stars with the highest possible spectral and 
spatial resolution. We consider here the case of NGC 4194, “The 
Medusa,” a galaxy with distorted morphology and luminous star 
formation that appears to be the result of an unequal merger 
between an elliptical and a smaller spiral. 

This is the fourth paper in a series on high-resolution 
spectroscopy of intense extragalactic star formation sources in 
the middle infrared (Beck et al. 2010, 2012, 2013). We use 
mid-infrared lines of metal ions to trace the kinematics and 
spatial distribution of ionized gas in Galactic H ii regions and in 
starburst galaxies (Achtermann & Lacy 1995; Zhu et al. 2008; 
Beck et al. 2013) because they are little affected by extinction 
and permit us to attain true spectral resolution, including thermal 
effects, much higher than is possible with any hydrogen line. The 
“Medusa galaxy,” NGC 4194 (a.k.a. Arp 160 and Mkn 201) is a 
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starburst galaxy at 39 Mpc distance (1" = 190 pc) with striking 
tidal features; it has been variously classed as Magellanic, 
Sm(pec), and BCG. The central kpc of NGC 4194 is a powerful 
infrared source with 60 /zm/100 /zm ~ 1, typical of starburst 
heating. Spitzer Infrared Spectrograph (IRS) spectra (from the 
Heritage Archive) show a spectrum typical of a starburst- 
dominated galaxy, with little if any contribution from an active 
galactic nucleus (AGN). NGC 4194 is believed to be a merger 
remnant, with the favored history being a small gas-rich spiral 
falling into an elliptical four times its mass (Manthey et al. 
2008). The central region of NGC 4194 is the site of intense 
star formation: Weistrop et al. (2004) and Hancock et al. (2006) 
find numerous bright optical and UV knots that they identify 
as very young globular cluster precursors. But the structure of 
this central starburst has not been probed. The observations are 
either too low in spatial resolution (e.g., the Two Micron Sky 
Survey) or too affected by the deep extinction (Hattori et al. 
2004) to determine the structure. Beswick et al. (2005) mapped 
the center of NGC 4194 at 21 cm with sub-arcsecond resolution 
and found two compact sources separated by only 0'/35, which 
they identify with the Galactic nucleus, and which Manthey 
et al. (2008) appear to tentatively accept as the nuclei of the 
progenitor galaxies, but the nature of these sources and their 
relation to the merger history have not been explored. 

We report here on radio continuum maps of NGC 4194 
with sub-arcsecond spatial resolution and on spectra of the 
12.8 /zm emission line of [Neil] with velocity resolution better 
than 5 km s _1 . The radio maps are at short cm wavelengths, 
which are sensitive to both the thermal emission of H ii regions 
and the non-thermal emission of supernova remnants. The fine- 
structure line of Ne + at 12.8 /zm is usually one of the strongest 
mid-infrared emission lines in Hu regions and is a preferred 
kinematic probe because of its low susceptibility to thermal 
broadening; in a T e = 7500 K Hu region the FWHM from 
thermal effects will be 4.1 km s -1 compared to 18 km s -1 
for H + . 
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Table 1 

Observations and Reductions 
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Date 

Instrument 

Wavelength 

Beam Size 

®las a 

Noise 

Program 

2013 Feb 2 

TEXES 

12.8 /zm 

1"4 

n.a. 

n.a. 

n.a. 

1987 Jul 6 

VLA 

2 cm 

O'.' 18 x (X'149 

376 

0.06 mJy bm -1 

AE47 

1987 Nov 21 

VLA 

6 cm 

1'.'4 x 1"13 

24" 

0.05 mJy bm -1 

AS286 

2000 Nov 20 

VLA 

6 cm 

0'.'27 x O'.' 17 

879 

0.08 mJy bm -1 

AN095 

” 

” 

” 

0'.'46 x o:'43 

879 

0.02 mJy bm -1 

” 

” 

” 

3.6 cm 

OC'18 x O'.' 104 

573 

0.023 mJy bm -1 

” 

” 

” 

” 

0'.'26 x 0'.'23 

573 

0.015 mJy bm -1 

” 

2001 Mar 19 

” 

” 

0?53 x 0:'53 

17" 

0.04 mJy bm -1 

” 


Note. a 0ias is the angular size of the largest structure to which the observations were sensitive. 


2. OBSERVATIONS 

2.1. TEXES [Neil] Data Cube 

NGC 4194 was observed at the NASA Infrared Telescope 
Facility on Mauna Kea, Hawaii, on the night of 2013 February 2, 
with the TEXES spectrometer (Lacy et al. 2002). TEXES on the 
IRTF has a seeing and diffraction-limited beam of ~174. The 
observations were taken in the high-resolution mode, with a 
plate scale 0736 x 0.9 km s -1 per pixel and spectral resolution 
~3 km s -1 . Ceres was the telluric comparison source and the 
wavelength scale was set from atmospheric lines. The slit was 
29 pixels long and l / . / 32 wide; it was set north-south and scanned 
across the galaxy in R.A. steps of 0 ".1 to generate data cubes 
with Nyquist spatial sampling. The total observation time was 
2.23 hr. 

2.2. Radio Maps 

NGC 4194 was observed at the NRAO Very Large Array 
(VLA) 7 in program AN095. A array observations at 6 and 3.6 cm 
(C and X bands) were obtained on 2000 November 17 and B 
array observations at X band on 2001 March 19. The data were 
calibrated and reduced with AIPS. Maps were produced with 
varied weightings and beam sizes; the observing and reduction 
parameters are given in Table 1 . 

2.3. Archival Data 
2.3.1. Radio 

We obtained archival data in the radio and infrared regimes 
from the VLA and from Spitzer. NGC 4194 was observed at 
2 cm (U band) in the A array for program AE47 and at 6 cm 
with the B array in program AS286. Note that the resolution 
of these B array observations is on the order of the diffraction- 
limited TEXES beam. The parameters of these maps are given 
in Table 1 . 

2.3.2. Infrared Spectra 

NGC 4194 was observed by Spitzer as part of the IRS 
Standard Spectra program and the results are in the Spitzer 
Heritage Archive. It is dominated by the low-excitation lines 
of [Nen] and [Sin] 18.7 pm. All the lines appearing in the 
9.9-19.6 pm SH and 18.7-37.2 pm LH modules are in the 
ratios expected of a starburst galaxy (Sturm et al. 2002); there 
is no spectral signature of an AGN. 


7 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 


3. Hu REGIONS IN THE CENTRAL Kpc OF NGC 4194 
3.1. Radio Continuum 

Optical and UV observations of NGC 4194 have detected 
multiple star-forming knots near the nucleus (Weistrop et al. 
2004; Hancock et al. 2006), but because of heavy dust obscu- 
ration these observations did not penetrate the nucleus itself. 
Published infrared and radio observations are either of low res- 
olution or at a wavelength not sensitive to the thermal emission 
of young star-forming regions. The radio observations we report 
here are the first high-resolution, obscuration-free images of star 
formation activity in the nucleus of NGC 4194. We first discuss 
the spatial distribution and structure of the nucleus, and then the 
nature of the radio sources. 

3.1.1. Spatial Distribution of Radio Emission 

Maps of NGC 4194 at 6 and 3.6 cm are shown in Figure 1, 
and a 2 cm map, converted from the B1950 to J2000 system, 
is in Figure 2. This is the first complete and extinction- 
free high resolution view of ionized gas structure near the 
nucleus in NGC 4194. The radio continuum in the central 
kpc of NGC 4194 has a complex structure which our highest 
resolution maps barely resolve. The central source 8 at a = 
12 h 14 m 9 s .68, 8 = +54 0 3E3578, appears elongated SE-NW at 
the lower resolutions. In the highest resolution 3.6 cm maps it 
breaks up into 3 sources in a line at p.a. ~ —40° . Another distinct 
source lies west of the center. South of the center is a region 
~2" in length, with three distinct sources in a background of 
extended clumpy emission. These seven distinct sources were 
fit with single Gaussians and their sizes, positions and fluxes are 
in Table 2. We discuss the nature of these sources in the next 
section. 

In addition to the compact sources there is extended emission 
over a 6" x 6" region. Figure 3 shows a naturally weighted 3.6 cm 
image plotted so as to emphasize the lower level emission. 
It is clearly not axisymmetric. Weak emission connects the 
central and western radio source with what Beswick et al. (2005) 
described as an “arm-like feature,” and the southern and central 
source are immersed in what could be a nuclear spiral, a nuclear 
bar, or part of a star-forming ring or disk. Non-axisymmetric 
structure has been seen in NGC 4194 and interpreted as a bar 
(Jog & Maybhate 2006), but that was at much larger (kpc) 


8 This source has been identified with the Galactic nucleus in earlier radio 
maps, and this has been assumed by other authors: Kaaret & Alonso-Herrero 
(2008), for example, attribute the X-ray source they find near here to the 
“nucleus.” However kinematic and other considerations, discussed in 
Section 4.3, cast doubt on that identification, and it will be called “the central 
source” in this paper. 
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Figure 1 . Uniform-weight maps at 6 (top) and 3.6 cm (bottom). For the 6 cm map, the B and A array observations were combined to create a data set with a O'/ 2 8 x O'/ 19 
beam, shown in the lower left, which was sensitive to structures as large 24". The lowest positive contour is 2.5cx = 0.095 mJy. The 3.6 cm map has a O'/ 18 x O'/I 
beam, shown, and the lowest positive contour is 2.5cx = 0.06 mJy. For both maps the levels are spaced logarithmically by factors of 2"/ 2 . 


radii and does not seem related to the nuclear emission. These 
structures are discussed in Section 4. 

3.1.2. Thermal Emission and the Spectral Index 

The radio continuum of starburst galaxies combines non- 
thermal emission from synchrotron sources created in previous 
generations of stars with thermal free-free emission generated in 
young H ii regions tracing the current episode of star formation 
(Condon 1992; Turner & Ho 1994). Non-thermal processes 


typically dominate the longest wavelengths (20 cm and longer) 
and thermal free-free the shorter cm and mm bands. The radio 
spectrum can be characterized by the spectral index a, such that 
S v = v a : for pure thermal free-free emission from H ii regions, 
a = —0.1, for optically thick emission a > —0.1 and can 
be as great as +2. For non-thermal synchrotron emission a is 
steeper, i.e., more negative, than the thermal value. Non-thermal 
emission of starburst galaxies is observed to have a in the range 
—0.5 and —1.2 between 21 and 2 cm. The integrated 21 cm 
flux of Condon et al. (1990) and the single dish 6 cm result of 


3 


The Astrophysical Journal, 787:85 (12pp), 2014 May 20 


Beck et al. 



Figure 2. Uniform-weight 2 cm map. The O'/ 18 x O'/ 15 beam is shown at lower left, the lowest positive contour is 0.15 mJy = 2.5a and the levels are spaced 
logarithmically by factors of 2 n / 2 . 


Table 2 

Source Parameters from 3.6 cm Maps a 


ID 

R.A. (2000) 

Decl. (2000) 

Major Axis 

Minor Axis 

Peak Intensity 

Total 

Central 1 

12 14 09.68 

54 31 35.73 

<y.'46 

0?22 

0.8 mJy bm -1 

1.2 mJy 

Central 2 

12 14 09.67 

54 31 35.74 

( y/21 

O'.' 14 

1.16 mJy bm -1 

2 mJy 

Central 3 

12 14 09.65 

54 31 35.86 

<y.'22 

O'.' 11 

1.9 mJy bm -1 

3.3 mJy 

W. Source 

12 14 09.326 

54 31 35.68 

a'23 

a' 14 

0.46 mJy bm -1 

0.94 mJy 

South 1 

12 14 09.59 

54 31 34.6 

<y.'5 

0'.'37 

0.3 mJy bm -1 

0.57 mJy 

South 2 

12 14 09.63 

54 31 34.27 

(y/51 

ff.'43 

0.46 mJy bm -1 

0.9 mJy 

South 3 

12 14 09.66 

54 31 34.21 

(y.'69 

0'.'513 

0.46 mJy bm -1 

1.24 mJy 


Note. a The nuclear and western cluster fits are from the uniform- weight map with a O'/ 18 x O'/I beam; the southern clusters are fit 
from the natural weighted map with a O'/ 26 x O'/ 2 3 beam. 


Becker et al. (1991) give the spectral index of NGC 4194 as 
—0.81 between 21 and 6 cm. 

The archival 2 cm observations offer a wider baseline view 
of the spectral index. We created a 2 cm image with the same 
beam size as the 6 cm and took the spectral index between 6 
and 2 cm, which is shown in Figure 4 along with the 2 cm map 
contours. The 6-2 cm spectral index shows that the central and 
southern sources have substantial non-thermal emission and the 
western is purely thermal. In fact, the spectrum of the western 
source may actually be rising at 2 cm, consistent with a thermal 
source of such high emission measure as to be partly optically 
thick at 6 cm. 

The radio continuum at 6 and 3.6 cm is a mixture of thermal 
free-free emission and non-thermal synchrotron emission, but 
the 2 cm map is mostly thermal. This is shown by comparing 
the fluxes to those at 2.7 mm, where the non-thermal emission is 
insignificant and dust emission is negligible. Aalto et al. (2010) 
measured 4.6 zb 1 mJy total flux at 2.7 mm over this entire region. 
With a — —0.1 this predicts total thermal fluxes 5.6 =b 1 mJy at 
2 cm and 6.3 zb 1 mJy at 6 cm. So 88% of the total 2 cm flux is 


thermal (but only 15% of the total 6 cm.) The 2 cm emission in 
Figure 2 is thus an excellent high-resolution representation of the 
extinction-free star formation traced by the free-free emission. 

The Lyman continuum rate to maintain this thermal flux is 
estimated from the prescription of Ho et al. (1990) and the 3 mm 
flux to be A^Lyc = 1 .0 x 10 54 (D / 39 Mpc) 2 s -1 for the three main 
star-forming regions of the central ~10". This corresponds to 
L 0 b = 3 x 10 10 L 0 , and a star formation rate of ~10 M 0 yr -1 
for a Kroupa 3 Myr starburst based on STARBURST99. 

The non-thermal emission in the compact sources of 
NGC 4194 is probably from recent supernovae. There must 
have been considerable supernova activity in the recent past: 
the 3.3 mJy of 3.6 cm emission which we calculate to be non- 
thermal requires 450 objects like Cas A. Beswick et al. (2005) 
found that at 21 cm, where the non-thermal component domi- 
nates, the emission appears with a MERLIN beam of O'/ 1 7 x O'.' 1 5 
to be a very compact double source. The two components have 
12.6 and 5.1 mJy and their positions and separation are con- 
sistent with Central Sources 1 and 3 in our 6 and 3.6 maps. 
There is also a weak 21 cm source north of the midpoint of 
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Figure 3. Low-level emission at 3.6 cm. The Ct/32 x (f/22 beam is shown. The lowest contour is 1.3cr (0.02 mJy bm ! ) and contour levels are 1, 1.4, 2, 2.8, 4, 6.25, 8, 
16, 32, and 64. 
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Figure 4. Spectral index between 6 and 2 cm in grayscale with the 2 cm intensity in contours. The beam is (f/26 x O'/ 2 3 and contour interval and levels as in Figure 3. 


the double source with peak flux density 0.7 mJy bm -1 ; it has 
no 3.6 cm counterpart. We cannot, with our present data, map 
the spectral index with sufficient spatial resolution to locate the 
non-thermal sources precisely. The 6-2 cm spectral index hints 


that the emission is more strongly non-thermal on the north- 
ern side of the central source, perhaps related to the weak third 
source in the 21 cm map, but that is not conclusive: the 2 cm 
data were observed in an epoch different from the 6 cm and the 
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Figure 5. Spatial distribution of the total line flux in the [Nen] data cube, in red. Contours are integer multiples of 10% of the peak flux of 9 x 10 -3 erg (s cm 2 sr) -1 . 
The apparent peak at R.A. 12 14 09.9 is spurious, from a flaw in the detector. Black contours show the 2 cm radio emission convolved to the [Nen] beam; contours 
are integral multiples of 1.68 mJy bin -1 and dashed contours are negative. 

(A color version of this figure is available in the online journal.) 


precession may not be perfectly accurate. We conclude that 
thermal and non-thermal emission are mixed in the central 
and southern sources of NGC 4194 at wavelengths 3.6 cm and 
longer. 

The purely thermal index of the small western source argues 
that it is a young Hu region that has not yet produced a 
significant number of supernovae. From the radio flux we 
deduce that the Lyman continuum flux in that source is A/flyc ~ 
2.08 x 10 53 y s _1 , equivalent to 2.08 x 10 4 07V stars and a 
star formation rate of ~2 M© yr _1 . If this is one star cluster it is 
one of the largest super star clusters yet seen; it is much brighter 
than the central source in NGC 5253 and almost as strong as 
the SW radio source in NGC 4102 (Beck et al. 2010) and the 
brightest source in M82 (Tsai et al. 2009). The extended nature 
of the western source argues that it is probably more than one 
cluster, however. 

In summary, the three main radio sources are intense regions 
of star formation, but at different evolutionary stages. The central 
and southern source are old enough that many of their stars have 
gone supernova. The presence of an X-ray source (presumably 
an X-ray binary) near the central radio source (Kaaret & Alonso- 
Herrero 2008) is another sign of stars in very late stages of 
evolution. The western source, in contrast, is too young to have 
formed a significant number of supernovae. 

3.2. [Nen] Spatial Structure Agrees with the Radio 

The spatial distribution of the [Ne ii] emission is found from 
the data cube by collapsing the cube along the velocity axis. 
The resulting map is shown in Figure 5, superimposed on a 
radio map convolved to the same beam size. In the [Ne n] map 


the western source is a discrete peak, while the southern source 
is a plateau extending the central source; the radio continuum 
at the same spatial resolution looks very similar. [Ne ii] is a star 
formation tracer showing the presence of young H ii regions, so 
its distribution here argues that all the nuclear sources, including 
those with non-thermal radio spectra, host some recent star 
formation. 

The total flux of [Nen] summed over the line extent and 
the region mapped is 1.6 x 10 -12 ergs -1 cm -2 , with absolute 
uncertainty ~30%. The IRS on Spitzer had two pointings of 
the SH module on this position and measured 4.7 and 5.0 Jy, 
which is 2 zb 0.06 x 10“ 12 erg s -1 cm -2 . There are considerable 
systemic uncertainties in converting Jy to flux for an unresolved 
line and in comparing observations of such different resolution 
as TEXES and IRS; we can say only that the results are 
consistent with TEXES having recovered the entire Spitzer flux. 

3.2.1. [N e + ] Abundance 

The strength of the [Ne ii] line and of the thermal radio con- 
tinuum are related because both depend on the emission mea- 
sure n 2 e l (cm -6 pc) If the neon abundance has the solar value 
of 8.3 x 10 -5 n(H) by number, and is all in Ne + (an overesti- 
mate, since Spitzer detected [Nem]), the collision strengths of 
Osterbrock & Ferland (2006) give 7 7 ([Nen])(erg s - 1 cm -2 ) = 
2.0 x 1 0 10 E 5 ©Hz (Jy)- This predicts a total 5 GHz flux of 
8 mJy for NGC 4194, compared to the 6 mJy deduced above for 
thermal radio emission. The sizable uncertainty in our absolute 
flux calibration permits us to say only that the flux is consistent 
with a solar or somewhat higher abundance, and even marginally 
consistent with the slightly less than solar metallicity cited by 
O’Halloran et al. (2006). 
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Figure 6. Radio continuum contours of the VLA 2 cm map are superimposed on the CO integrated intensity and CO velocity field from Aalto & Huttemeister (2000). 
The estimated accuracy of the overlay is ^(f/2. The velocity contours are 2400 (blue) to 2640 (red) km s -1 with 24 km s -1 spacing. The CO contours are the total 
integrated intensity of CO emission with levels (0, 1.5, 3.5,...) Jy bm -1 km s -1 . Aalto & Huttemeister (2000)’s peaks A, B, C, and D go clockwise from the southeast 
and E is the northeast feature. 

(A color version of this figure is available in the online journal.) 


3.3. Ionized Gas Does Not Coincide with the Molecular Gas 

Aalto & Huttemeister (2000) used the OVRO interferometer 
to map the 12 CO(1-0) emission in the center of NGC 4194 
with beam sizes l'/7-2'.'5 and velocity resolution of 20 km s _1 . 
They found that in the region mapped in radio and [Neil] 
the molecular gas comprises 5 distinct peaks and an extended 
component. We precessed their figures to the J2000 system used 
in the radio continuum observations and show the 2 cm map 
overlaid on the CO in Figure 6. Since both the radio and CO maps 
are interferometer maps, the uncertainty in registration should be 
negligible. It is clear from the figure that the spatial distribution 
of molecular gas differs materially from the ionized. First, the 
brightest radio and [Nen] source is not perfectly coincident 


with any of the five CO peaks, but sits on the northwest side 
of CO peak B. The offset is ~0 / . / 6, much larger than the 071 
conservatively estimated for the possible error in the precession. 
Second, the southern clump of radio sources agrees in position 
with the molecular peak A, but not in relative intensity: the 
radio and [Ne ii] flux on peak A is much weaker than on peak 
B, which has the brightest radio emission. The CO/2 cm ratio 
on A is higher than on B by at least a factor of five. Third, there 
is no ionized gas detected on peak D even though it is almost 
(80%) as strong in CO as peak C, which is coincident with the 
western radio source; the CO/2 cm ratio in D must be at least 
10 times higher than in C. 

We see that the greatest concentration of molecular mass is 
neither the radio peak nor the geometrical center of the starburst 
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activity, and that star formation activity is concentrated on one 
side of the molecular gas distribution. In other words, the ionized 
gas does not trace the mass distribution on scales of ~0 f /5 or 
^80-100 pc. The mass of young stars in the Hu regions is not 
significant on these scales: the total stellar mass in the clusters 
is less than a few times 10 7 M 0 , a small fraction of the total 
molecular mass of «2 x 10 9 Mq and of the «10 8 M 0 in each 
cloud. 

This raises a new question: where and what is the Galactic 
nucleus?. We cannot assume that the nucleus is identical with 
one of the radio peaks — not even the central source. We return 
to this in the Conclusions. 

4. KINEMATICS 

The kinematics of molecular gas in NGC 4194 were studied 
by Aalto & Huttemeister (2000), who observed CO emission 
with a 2'.'0 x 1'.'7 beam and 20 km s -1 spectral resolution, and 
by Aalto et al. (2010) who measured 13 CO with a 4'.'56 x S'.' 9 8 
beam. The atomic gas was observed by Beswick et al. (2005), 
who had 21 km s -1 resolution. The velocity field in NGC 4194 is 
extremely complex, reflecting the merger that created the galaxy. 
Aalto & Huttemeister (2000) and Beswick et al. (2005) agree 
that in the inner ~0.65 kpc (3'.'5) there is a shallow N-S velocity 
gradient of ~60 km s -1 arcsec -1 , consistent with rotation, that 
there are significant non-circular motions at larger distances, 
and that the nuclear region contains some velocity components 
not due to rotation. 

4.1. Kinematics of Ionized Gas: [Nen] Distribution 

The TEXES observations are the first to measure ionized gas 
kinematics in NGC 4194. The spatial distribution of the emission 
is shown as a two-dimensional grid of spectra in Figure 7(a). The 
[Ne n] emission in the central 5" x 5" covers the velocity range 
2375-2650 km s -1 . This is close to the full range over which 
Beswick et al. (2005) see Hi absorption in the central 3". The 
overall distribution of the peak velocities also agrees roughly 
with the H i result: the highest velocities observed are around 
2600 km s -1 and are in the southern extension, the lowest are 
around 2400 km s -1 and are at the north end of the source, and 
the northwest is at intermediate velocities of ~2500 km s -1 . 

Besides the similarities there are substantial differences 
between the ionized gas kinematics, as shown on the spectral 
grid, and the atomic and molecular findings. While Ne + is 
the commonest ionization state of neon in a high-metal star 
formation region like NGC 4194 and is probably widely 
distributed in the galaxy, the [Ne ii] line is collisionally excited 
and is most readily seen in Hu regions of moderate to high 
density (i.e., >10 3 cm -3 ). So the [Nen] is concentrated in the 
compact and intense radio sources to which the high-resolution 
radio maps are sensitive. The spatially isolated 2500 km s -1 
feature may be straightforwardly associated with the western 
thermal source. We assign the 2440 km s -1 feature to the 
nuclear source and the 2600 km s -1 to the southern emission 
complex. The spatial resolution is not high enough to separate 
the central and southern sources completely and they overlap 
in some positions, creating very complex line profiles (as for 
example, in Figure 7(b)). 

4.2. [Nen] Line Profiles 

What can the [Nen] spectra tell us about the ionized gas 
motion in this galaxy? The profile of a fine- structure emission 
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line can be formed and influenced by thermal velocity disper- 
sion, gravitational turbulence in a virialized system (such as a 
bound nebula), and non-turbulent bulk motions. The obvious 
bulk motion in NGC 4194 is rotation creating a velocity gradi- 
ent. Does this contribute to the line profiles, or can we treat the 
sources as kinematically isolated? We explored this question 
by creating data cubes from the spatial structure of the radio 
maps and a given velocity gradient. For a reasonable gradient 
like the ~60 km s -1 arcsec -1 of the CO, the gradient does not 
affect the line widths of the individual sources. The sources are 
simply too small to include a significant range of the gradient 
velocities. This justifies analyzing the line profiles in terms of 
internal, rather than large-scale Galactic, motions. 

Thermal velocity dispersion is not an important factor in 
these line profiles; because of the high mass of neon an Hu 
region with T e = 7500 K has FWHM from thermal effects 
of only 4.1 km s -1 . The FWHM increases only as the square 
root of T e . Gravitational turbulence in a virialized system will 
create a Gaussian line whose FWHM can be found from 
o = FWHM/2.35 = {GM /3R) 1 ^ 2 where R is the radius of 
the system. But since we do not resolve the individual emitting 
regions in the [Ne ii] sources we do not know R. Further, the 
radio maps show that the central and southern [Nen] source 
each include several distinct radio sources, presumably super 
star clusters, whose sum creates the observed lines. So while the 
gravitational turbulence is probably responsible for most of the 
width of the line profile, it cannot be satisfactorily calculated 
from these data for the southern and central sources. For the 
spatially isolated and apparently simple western source, we 
show the line summed over a 3 x 3 pixel box in Figure 8 along 
with the best single Gaussian fit. The FWHM is 104 km s -1 
and the line center is 2480 km s -1 ; x 2 /d.o.f. is 1-03 and 
the residuals do not show structure. The virial velocity relation 
formally gives M Q /R( pc) ~ 1.4 x 10 6 M 0 pc -1 , while the total 
stellar mass estimated from the Ai yc and a Kroupa initial mass 
function (IMF) extending to 0.1 M 0 is only 8 zb 3 x 10 5 Mq. 
The source is unresolved with the O'.' 18 (~35pc) beam of the 
highest resolution radio map. It is likely that even this apparently 
simple source contains multiple clusters at offset velocities, 
which make the line appear wider than the true turbulent 
FWHM. Fine profiles in the central and southern sources have 
FWHM from 85-104 km s -1 , which formally gives M/R of 
0.8-1. 4 x 10 6 M o pc -1 , but since these sources clearly (from 
the radio maps) contain multiple sub-sources both the FWHM 
and the mass concentrations are only upper limits. 

4.3. [Nen] Kinematics on Large Scales: 
Position-Velocity Diagrams 

To examine what [Nen] can tell us about the large-scale 
kinematics of the nuclear region, we created a data cube 
binned by 4 spectral pixels and 2 pixels along the slit, to 
create almost square O'.' 7 x O'.' 7 x 3.74 km s -1 pixels. From 
this we can see the ionized gas kinematics with minimal 
spatial averaging. Figure 9(a) shows the line extent as a 
function of position along a N-S cut through the central source; 
Figure 9(b) shows an E-W cut, also through the center. The 
N-S position-velocity diagram has an obvious velocity gradient 
of several tens of km s -1 arcsec -1 with blue north and red 
south, very similar to the CO (Aalto & Huttemeister 2000) and 
Hi (Beswick et al. 2005) results. The ionized gas gradient is 
200 km s -1 in 4'.' 2, for a spatial gradient (corrected for the 40° 
inclination angle) of 390 km s -1 kpc -1 . The total velocity extent 
is similar to the CO results but the gradient is slightly higher 
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Figure 7. Top: the observed [Nen] spectra in each pixel. The intensity units are erg (s cm 2 sr) -1 . Bottom: sample pixels from the above figure which are enlarged to 
show the overlapping velocity components. 


than the 320 km s -1 kpc -1 of the atomic and molecular gas; this 
may be because of the very different beam sizes. 

The [Ne n] observations average over a very patchy distribu- 
tion of ionized gas. We checked that this does not distort the 
position-velocity diagram by comparing it to that of the sim- 
ulated data cube based on the radio continuum maps; they are 
consistent. The gradient is not due to averaging over sources at 
offset velocities: the velocity offsets between the [Ne n] sources 
result from their positions along an incompletely sampled 
velocity gradient. 

In the E-W position-velocity diagram the central and the 
western source appear spatially separated and with almost 


identical velocity ranges of 2400-2550 km s _1 . There is nothing 
to suggest an E-W gradient; we find from the simulated cube 
that if there were a gradient as large as 20 km s -1 arcsec -1 , 
it would appear in the position-velocity diagram. There is a 
separate red feature at 2580 km s -1 , most clearly in the central 
source; we believe that to be part of the southern source which 
was included in the same beam. 

4.4. A Galaxy Core in Solid-body Rotation 

In the previous sections we showed that the none of the 
H ii regions measured in the [Ne ii] and radio are likely to be 
a Galactic nucleus. But even if the mass distribution is not 
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Figure 8. Line profile of the entire western source, with the best single Gaussian fit and the residuals. The fit center is 2480 km s 1 , FWHM 103.9, and x 2 /dof = 1 .03. 
(A color version of this figure is available in the online journal.) 




Figure 9. Position-velocity diagrams through the [Nen] sources. For these diagrams the data was binned by 4 pixels in velocity and by 2 pixels along the slit; the 
final scale is O'/ 7 x O'/ 7 x 3.7 km s _1 . Top is an N-S cut through the nucleus and N is up; bottom an EW cut through the brightest source with E up; in both, blue is to 
the left. 
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apparent in the CO, radio, and [Neil] maps, its gravitational 
effects cannot be hidden. We now find that there is no kinematic 
sign of a nuclear mass concentration either. In a spiral galaxy 
there is usually a very broad line on the nucleus and a 
velocity jump across it; there is no such pattern in NGC 4194. 
Specifically, the central source, the apparent candidate for the 
nucleus, has the narrowest line profile of all the [Ne n] sources. 
The only kinematic hint at a mass concentration is in the CO 
velocity field; the isovelocity contours, which are evenly spaced 
across the main radio and CO peaks, kink and take on a closer 
spacing around a = 12 h 14 m 09 s .5, 8 = +54°3E36". But that 
location is between the central and western radio sources and 
there is no emission peak, at any wavelength, there. Instead, 
the kinematics are dominated by a smooth gradient N-S in 
both molecular and ionized gas. Smooth gradients imply solid- 
body rotation. In disk galaxies solid-body rotation is usually 
associated with rings, and elliptical galaxies, especially those 
that are post-merger, may host cores several hundred pc in 
size, kinematically distinct from the rest of the galaxy, with 
this rotation pattern (de Zeeuw & Franx 1991). 

How can the galaxy have the kinematic signature of solid- 
body rotation while the ionized gas appears to be a two-armed 
spiral? If we critically examine the low-level radio emission 
maps in Figure 3, it is clear that the appearance of “spiral 
structure” depends very strongly on the southeast elongation 
of the southern source. The sub- sources fall on a line southeast- 
northwest and the emission extends further in a spur in the same 
direction. The whole southern complex is almost perpendicular 
to the emission connecting the southern and central sources, 
so that together with the elongated central source the region 
suggests a “twin peaks” star formation structure — but again, 
that contradicts the kinematics, which are not at all like those of 
the barred spirals where “twin peaks” form. Finally, the “arm- 
like” arc of radio emission connecting the central and western 
sources lies northwest of the main CO emission in the region of 
highest radio-to-CO ratio; it may show a front of star formation 
proceeding into the clouds, rather than a true spiral arm. So the 
spiral structure which appears to first glance is less certain on 
closer inspection. 

Let us now work from the opposite direction: how does the 
ionized gas distribution fit the kinematic results? The central 
and southern sources are clusters that have formed in a rotating 
spheroidal core with southern end receding. We cannot rule out 
that the rotating region is a very thick ring or disk instead of a 
spheroid; if this is the case the clusters have formed throughout 
the thickness of the ring. That the western source has almost the 
same velocity as the central may agree with the finding of Aalto 
& Huttemeister (2000) that the velocities of clouds C and D are 
consistently lower than if they were in solid-body rotation with 
A and B. The western source, like its host cloud, may not have 
yet settled into equilibrium. 

5. CONCLUSIONS 

We have presented high-resolution radio continuum maps at 
6, 3.6, and 2 cm, and a high spectral-resolution velocity-position 
cube in the [Nen] 12.8 /zm emission line, of the central 
1.3 kpc of NGC 4194. These maps are the first high-resolution, 
extinction-free, measurements of ionized gas in the center of 
this galaxy and trace the star formation regions with resolution 
up to 0? 18 (~35 pc). We find that: 

1 . NGC 4194 hosts a nuclear starburst that has formed multiple 

compact sources, apparently groups of embedded super star 
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clusters, within a 3" radius. The low-level radio emission 
suggests a “nuclear spiral” but that is not confirmed by the 
kinematics. 

2. Based on the thermal radio emission at 2 cm, which 
agrees well with previous observations at 2.7 mm, we find 
for the central ~10" region a thermal free-free flux of 
^cm™' = 5.6 zb 1 mJy, which implies a Lyman continuum 
rate of N^ yc = 1.0 x 10 54 (D/39Mpc) 1 2 s -1 . For a Kroupa 
IMF and a 3 Myr starburst, this implies a star formation 
rate of ~10 M© yr -1 , and a luminosity in massive young 
stars of Lob = 3 x 10 10 L 0 , about a third of the total 
infrared luminosity of Lir = 8.5 x 10 10 L o . The star 
formation in this galaxy is dispersed widely over the 
inner 10" (1.9 kpc) region. Our radio observations are not 
sensitive to regions with 2 cm fluxes less than ~0. 1 mJy, or 
Ar y C < 2 x 10 52 (D/39Mpc) 2 s -1 , so it is possible that the 
remaining infrared luminosity arises in smaller, and widely 
dispersed, star- forming regions within the system. 

3. Of the compact emission sources that are bright enough to 
analyze, one has a pure thermal or slightly rising spectrum, 
typical of a very young embedded H ii region. The others 
are older and have a mix of thermal and significant non- 
thermal emission. The [Ne ii] line velocity dispersions are 
consistent with the gravitational effects of many dense star 
clusters. 

4. The ionized gas has a smooth N-S velocity gradient of 
390 km s -1 kpc -1 across the observed region, consistent 
with a core in solid-body rotation. 

5. None of the radio sources can be convincingly identified 
with a Galactic nucleus. The nuclei of the original merger 
partners have not been identified. 

6. The CO velocity field shows that a mass concentration is 
present at a = 12 h 14 m 09 s .5, 8 = +54°3E36". This location 
is in between and roughly equidistant from the three star 
formation regions, and there is no star formation detected 
there. 

How can the last two items on this list be reconciled with the 
scenario that the Medusa is a minor merger? Each of the merger 
partners presumably hosted a nucleus of mass 10 6 -10 8 M© that 
would be detected via their influence on the kinematics. It is 
possible, and would be consistent with the simulations, to have 
both original nuclei now included in the mass concentration 
detected in the isovelocity contours. But in this case the nuclei 
must be devoid of star formation and ionized gas, or they would 
be detected in the radio and [Nen]. The gas and stars in the 
center of NGC 4194 have not come to equilibrium and the 
starburst may continue in regions now quiescent. Measurements 
of ionised and molecular gas with higher spatial resolution, and 
simulations exploring the central kpc, could perhaps find more 
clues to the behavior of this remarkable galaxy. 

TEXES observations at the IRTF were supported by NSF 
AST-0607312 and by AST-0708074 to Matt Richter. This 
research has made use of the NASA and IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion 
Laboratory, Caltech, under contract with NASA. We thank an 
anonymous referee for careful and thoughtful comments. 
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